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^ ■ Abstract 

■ We present a complete study of AS = 2 and AB = 2 processes in a warped extra 

■ dimensional model with a custodial protection of Zb]J)L, including ek, AMk, 
AMs, AMd, ^gL, AFg, Acp{Bd i)Ks) and Acp{Bs ipcp). These processes 

■ are affected by tree level contributions from Kaluza-Klein gluons, the heavy KK 



photon, new heavy electroweak gauge bosons Zh and Z', and in principle by 
tree level Z contributions. We confirm recent findings that the fully anarchic 
^ approach where all the hierarchies in quark masses and weak mixing angles are 

■ geometrically explained seems implausible and we confirm that the KK mass scale 

Mkk generically has to be at least ~ 20TeV to satisfy the ek constraint. We 
point out, however, that there exist regions in parameter space with only modest 
fine-tuning in the 5D Yukawa couplings which satisfy all existing AF = 2 and 
electroweak precision constraints for scales Mkk — 3TeV in reach of the LHC. 
Simultaneously we find that Acp{Bs — > ^/;(/>) and A^j^ can be much larger than in 
the SM as indicated by recent results from CDF and D0 data. We point out that 
for Bfi s physics AF = 2 observables the complex {Zh, Z') can compete with KK 
gluons, while the tree level Z and KK photon contributions are very small. In 
particular we point out that the Zd^^d?j^ couplings are protected by the custodial 
symmetry. As a by-product we show the relation of the RS flavour model to the 
Froggatt-Nielsen mechanism and we provide analytic formulae for the effective 
flavour mixing matrices in terms of the fundamental 5D parameters. 



1 Introduction 



The Standard Model of particle physics is in spectacular agreement with everything we 
know about interactions of elementary particles. Yet it requires large hierarchies to be 
put in by hand. There is no explanation for the hierarchy between the electroweak (EW) 
scale and the Planck scale and for the observed hierarchical pattern of fermion masses 
and mixings. 

Among the most ambitious proposals to explain these hierarchies are models with a 
warped extra spatial dimension first proposed by Randall and Sundrum (RS) [1] where 
the SM fields, except the Higgs boson, are allowed to propagate in the bulk [2-4]. These 
models provide a geometrical explanation of the hierarchy between the Planck scale 
and the EW scale and one can naturally generate the hierarchies in the fermion mass 
spectrum and mixing angles [2, 4] while simultaneously suppressing fiavour changing 
neutral current (FCNC) interactions [5,6]. Recently realistic models of EW symmetry 
breaking (EWSB) have been constructed [7-12] and one can even achieve gauge coupling 
unification [13,14]. 

In this work we discuss the fiavour structure of models based on the bulk gauge group 
Gbuik = SU{3)c X SU{2)l X SU{2)r x U{1)x x Plr. (1.1) 

The SM fermions are embedded in representations of Gbuik, so that there is a protection 
of the T parameter [7,8] and the coupling ZbibL [15]. This allows KK masses of order 
Mkk ^ (2 - 3) TeV which are in the reach of the LHC [10, 11, 16-18]. 

The goal of the present paper is to analyse the well measured FCNC processes related 
to particle- antiparticle mixings — and B^^^ — B^g. The off-diagonal mixing am- 
plitudes (i = K, d, s) receive dangerous tree level contributions from Kaluza-Klein 
(KK) gluon and EW gauge boson exchanges [6,19]. We would like to know whether this 
model can be made consistent with simultaneous constraints from Ek, AM/^ , AM^, AM^ 
and the mixing induced CP-asymmetry S^Ks for KK scales as low as Mkk — (2 — 3) TeV. 

A recent study [20], which applied model-independent results of the UTfit group [21] 
to RS-type models, concluded that the measured value of Ek implies that the lightest 
KK gluon mode has to be generically heavier than ~ 21 TeV, if the hierarchy of fermion 
masses and weak mixings is solely due to geometry and the 5D Yukawa couplings are 
anarchic and of 0{1). KK particles that heavy undermine the basic motivation for RS 
models. We would like to investigate if the KK scale could be lowered to be in reach 
of the LHC by allowing for a modest hierarchy and some tuning in the fundamental 5D 
Yukawa couplings. 

The dominant fiavour constraint comes from the CP-violating contribution to chi- 
rality fiip operators Qlr which are very strongly suppressed in the SM, but present in 
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RS models. The lower bound on the KK gluon mass obtained in [20] originates from 
the excessive contribution of Qlr to Ek- One of our strategies will be to find regions of 
parameter space consistent with EW precision observables [10,11,16-18] for which Qlr 
is sufficiently suppressed even if the 5D Yukawa couplings are mostly anarchic. 

Several alternative models have been proposed to deal with the flavour problem of RS. 
All depart from the fully anarchic set-up by incorporating some sort of flavour symmetry. 
One approach is to protect the model from all tree level FCNCs by incorporating a full 
5D GIM mechanism [22]. The bulk respects a U{3Y flavour symmetry and all flavour 
mixing is generated by kinetic terms on the UV brane. Although this model is safe, 
since its effective theory is minimal flavour violating (MFV) [23-27], it leaves the origin 
of the large hierarchies in the flavour sector unanswered. More recent proposals therefore 
seek to suppress dangerous FCNCs and simultaneously try to explain the hierarchical 
structure of the flavour sector. One of them is the so called "5D MFV" model [28]. 
Here one postulates that the only sources of flavour breaking are two anarchic Yukawa 
spurious. The low-energy limit is not MFV, and the additional assumption, that brane 
and bulk terms in the down sector are effectively aligned, is needed to suppress dangerous 
FCNCs. Recently, an economical model has been proposed [29] in which one assumes a 
f/(3) flavour symmetry for the 5D fields containing the right handed down quarks. This 
global symmetry forces the couplings of the right handed down quarks to the vector 
KK modes to be degenerate. Dangerous contributions to Qlr are only generated by 
suppressed mass insertions on the IR brane where the symmetry is necessarily broken 
(see [30] for a discussion of possibly problematic fermionic brane kinetic terms). Another 
recent approach [30] presents a simple model where the key ingredient are two horizontal 
U{1) symmetries. The SM fields are embedded into the 5D fields motivated by protecting 
ZhiJ)L- The horizontal U{1) symmetries force an alignment of bulk masses and down 
Yukawas which strongly suppresses FCNCs in the down sector. FCNCs in the up sector, 
however, can be close to experimental limits. In the present paper, however, we will 
study the original version of the model. 

As there have been other analyses of particle-antiparticle mixing in the RS model in 
the past [6, 19,31-34], most recently in [20], it is mandatory for us to state what is new 
in our paper: 

• First of all we perform a simultaneous analysis of the most interesting AF = 2 
observables in the K and Bd^s meson systems in conjunction with ek- In [20] only 
one Wilson coefficient at a time has been considered. This will give us a global 
picture of correlations between various observables. Such an analysis has not been 
performed in the literature so far. 

• Similarly we perform the full renormalisation group analysis at the NLO level. 
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including not only the two Qlr operators in our analysis, but also Qll and Qrr. 
We would like to emphasise that this is essential since the operator Qll, although 
subleading in ek, turns out to be as important as Qlr in Bd,s physics observables. 
On the other hand Qrr is subdominant in all processes considered in this paper. 

• In addition to tree level KK gluon exchanges considered in [20] we present for 
the first time the formulae for the EW tree level contributions of Z, Z', Zh and 
the KK photon A^^^ to the AF = 2 Wilson coefficients of the operators involved. 
Quite unexpectedly we find that the Z' and Zh contributions, while subleading 
with respect to KK gluon contributions in the case of Ek and AM^, can compete 
with the latter in the case of Bs^d physics observables. The contributions of the 
KK photon turn out to be small. 

• We point out and demonstrate explicitly that in the model in question tree level 
fiavour violating Z couplings to left-handed down-type quarks are strongly sup- 
pressed by the Plr symmetry up to small symmetry breaking effects due to the 
UV boundary conditions. This suppression mechanism works not only for the 
KK gauge boson contribution, but also for the KK fermion contribution to the Z 
coupling, as the fermion representations are symmetric under Plr- 

• We show that it is possible to simultaneously fit the SM quark masses and CKM 
parameters within their experimental 2a ranges and obtain agreement with all 
available constraints on AF = 2 observables. 

• We present a new useful parameterisation of the 5D Yukawa coupling matrices, 
taking into account only physical parameters. 

• As a by-product we analyse the connection of RS models to the Froggatt-Nielsen 
mechanism [35] and provide analytic formulae for the effective flavour mixing ma- 
trices in terms of the fundamental 5D parameters. 

Our paper is organised as follows. In Section 2 we summarise briefly the main ingre- 
dients of the model in question. Readers familiar with RS models may skip this section 
and start directly with Section 3, where we analyse the connection of flavour in RS mod- 
els with the Froggatt-Nielsen mechanism and derive explicit formulae for the effective 
flavour mixing matrices Ul,r and T>l,r in terms of the fundamental 5D parameters. Then 
in Section 4 we derive the effective Hamiltonians for — K^, — B^ and B^ — B^ 
mixings originating from tree-level KK gluon exchange and we calculate most interest- 
ing observables such as the CP-violating parameter ek, the mass differences AMk, AMd 
and AMg, the CP-asymmetries A%^^ (g = ci, s), Ac^{Bd — > ipKs) and Ac^{Bs — > ipcj)) 
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and the width difference ATq. We also give formulae for EW tree level contributions and 
estimate their size. Interestingly in the case of B^^s physics observables some of these 
contributions can compete with the KK gluon contributions. We also demonstrate that 
tree level flavour violating Z couplings to left-handed down-type quarks are strongly 
suppressed in the model in question, due to the custodial protection mechanism. This 
finding has important implications not only for AF = 2 processes, but also for AF = 1 
rare decays. In Section 5 we outline our strategy for the numerical analysis, presenting 
in particular a useful parameterisation for the 5D Yukawa couplings. In Section 6 a 
detailed numerical analysis of particle-antiparticle mixing observables is presented. We 
summarise our results in Section 7. 



2 The Model 

2.1 Geometric Setup 

The class of models we are considering is based on the Randall-Sundrum (RS) geometric 
background, i. e. we consider a 5D spacetime, where the extra dimension is compactified 
to the interval < y < L, with a warped metric given by [1] 

ds^ = e-^^^ri^^dx^'dx'' -dy^. (2.1) 

Here the curvature scale k is assumed to be ~ 0{M-p\). Due to the exponential warp 
factor e~^y , the effective energy scales depend on the position y along the extra dimension. 
In order to obtain a solution to the gauge hierarchy problem, we set e^^ = 10^^ and treat 

/ = ke-''^ ~ C(TeV) (2.2) 

as the only free parameter coming from space-time geometry. 



2.2 KK Gluons 



The main actors in AF = 2 processes in RS are at first sight KK gluons originating 
from the bulk SU{3)c symmetry, and in particular the first KK excitation. Therefore let 
us restrict ourselves to this mode, to be simply called KK gluon, in what follows. The 
profile of the KK gluon along the extra dimension is given by [2] 
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where Ji(x) and Yi{x) are the Bessel functions of first and second kind, 6 ~ is deter- 
mined by the boundary conditions at ?/ = 0, L and is a normalisation constant. The 
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KK mass Mkk can be numerically determined to be [36] 

Mkk^2.45/, (2.4) 

with / being defined in (2.2). We would like to caution the reader that a different 
notation has been used in [45]: Their Mkk corresponds to our /, so that in spite of 
comparable Mkk the masses of the first KK gauge bosons in that paper are larger than 
in our analysis. 

It is crucial to note that the KK gluon bulk profile is not fiat along the extra dimen- 
sion, but due to the factor e''^ strongly localised towards the IR brane. This localisation 
will give rise to flavour non-universal couplings of the KK gluon and eventually to tree 
level FCNC transitions, as we will discuss later on. 

2.3 Electroweak Gauge Sector 

The other actors in our analysis are the neutral EW gauge bosons Z, Zh, Z' and the KK 
photon [37], originating from the enlarged gauge group 

SU{2)l X SU{2)n x U{l)x x Plr . (2.5) 

We give here for the first time their contributions to AF = 2 processes using the couplings 
worked out in [37]. While subleading in the case of Ek and AM^ they turn out to play 
an important role in the case of B^^s physics observables. 

2.4 Bulk Fermionic Zero Modes 

Bulk fermions in RS models offer a natural explanation of the observed hierarchies in 
fermion masses and mixings [2,4,5] and provide at the same time a powerful suppression 
mechanism for FCNC interactions, the so-called RS-GIM mechanism [6]. 

The bulk profile of a fermionic zero mode depends strongly on its bulk mass parameter 
Cij,. In case of a left-handed zero mode \E'^^ it is given by [2,4] 

with respect to the warped metric. Therefore, for Cij, > 1/2 the fermion ^f^^"* is localised 
towards the UV brane and exponentially suppressed on the IR brane, while for cj, < 1/2 
it is localised towards the IR brane. The bulk profile for a right-handed zero mode ^'^^ 
can be obtained from 

fS\y,c^) = lf\y,-c^), (2.7) 
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so that its localisation depends on whether < —1/2 or c^. > —1/2. Note that the left- 
and right-handed zero modes present in the spectrum necessarily belong to different 5D 
multiplets, so that generally c*^ 7^ c,^^. 

In order to reproduce the SM quark content in the low energy limit, three left-handed 
zero mode SU{2)l doublets with bulk mass parameters Cq and three right-handed 
zero mode up- and down-type SU (2) ^ singlets u-Jj and d-Jj with bulk mass parameters c;^ ^, 
respectively, are required. As the KK fermions do not enter directly tree level AF = 2 
processes, we do not specify the fermion representations here. They are discussed in detail 
in [37]. We stress however that in order to preserve the custodial Zd\d-'j^ protection, Plr 
symmetric fermion representations are required. 

The coupling of a zero mode fermion '^'f^n to the KK gluon G^/^" in the flavour 
eigenbasis is given by 

„5D rL 2 

^tRG^^^^'^tn ■■ -^Yt^§j-,ldye'y[ff},iy,c^)\g{y), (2.8) 

with gf-' being the 5D SU{3)c gauge coupling constant and the SU{3)c generators. 
The 4D QCD coupling constant gf^ is then given by 

1 5D 

^f = ^^, (2.9) 

where in the absence of brane kinetic terms puv = 1- 

Note that flavour universality is strongly violated due to the dependence of the overlap 
integral on the bulk mass parameter Cij,. 



2.5 Brane Kinetic Terms 

One should keep in mind that localised brane kinetic terms for the gluon can change 
the relation between the bulk gauge coupling g^^ and the 4D QCD gauge coupling. The 
matching relation at a given scale contains both a bulk term and contributions from 
IR and UV brane kinetic terms. The UV brane terms consist of a possible (positive) 
bare kinetic term and an always present negative term [13,38-44], which is due to the 
asymptotically free QCD running from the Planck scale to the TeV scale, see [20] for an 
extended discussion. The running of the IR brane kinetic term is negligible and we will 
therefore focus on the UV brane localised kinetic terms. One possibility is that there are 
no bare UV brane kinetic terms at the Planck scale and that the negative contribution 
from the running reduces the bulk gauge coupling to gf^\/k ~ 3, corresponding to 
Puv ~ 1/2. Another possibility is that large brane kinetic terms at the Planck scale are 
present which would render the bulk strongly coupled, g^^y/k ~ Att. Usually one assumes 
an intermediate scenario where the bare kinetic terms are of exactly the same size as 
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the one induced by the one-loop running. This cancellation is assumed in the trivial 
matching relation (2.9) for puv = 1 and results in gf^y/k ~ 6 for the QCD bulk gauge 
coupling. In our analysis we have set puv = 1 which in [20] resulted in Mkk > 21 TeV 
in order to be consistent with the experimental value of Ek- 

Similarly, also in the electroweak sector, brane kinetic terms can be present and thus 
alter the simple tree level matching condition g^^ = g^^/y/L. Here the situation is 
additionally complicated by the fact that on the UV brane the gauge group is broken to 
SU{2)l X f/(l)y. Therefore different UV brane kinetic terms can be present for SU(2)l 
and U{1)y, so that in general p^y^'^^^ 7^ Puv^^ 7^ 1 ^^'^ a-lso different from puv iii the 
strong sector. Consequently the relative size of strong and electroweak contributions 
to AF = 2 processes depends on the size of possible brane kinetic terms. In order 
to allow for an easy comparison of the KK electroweak gauge boson effects with the 
KK gluon effects, we will, as in the former case, assume also here the intermediate 
scenario Puv'^^''^ ~ ^'uv^^ ~ order to keep the analytic expressions in Section 4.3 

simple, we omit Puv*'^''^'Puv''^ ^he formulae. A generalisation to include these terms 
is straightforward. 



2.6 Higgs Field and Yukawa Couplings 

The present analysis does not require the specification of the exact EWSB mechanism. 
Instead we will simply assume the presence of a Higgs doublet H{x^) residing on the IR 
brane. Once its neutral component develops a VEV v ~ 246 GeV, EWSB is achieved. 

The effective 4D Yukawa couplings, relevant for the SM fermion masses and mixings, 
are then given by 

y:^'' = Kf ^^f?\y = L, c},)fP{y = L, <j ^ a:;^'^ > (2-10) 

where A"''^ are the fundamental 5D Yukawa coupling matrices. In order to preserve 
perturbativity of the model, we require, as usual, \X^f\ < 3. Here and in the following, 
we work in the special basis in which the bulk mass matrices are taken to be real and 
diagonal. Such a basis can always be reached by appropriate unitary transformations in 
the Qi, Ui and di sectors. 

Due to the exponential dependence of F"''^ on the bulk mass parameters CQ^u,d, the 
strong hierarchies of quark masses and mixings can be traced back to 0{1) bulk masses 
and anarchic 5D Yukawa couplings A"'*^. We will elaborate more on this issue in the next 
section. 

The transformation from the quark flavour eigenbasis ul,r, dL,R to the mass eigenbasis 
ul,Ri dL,R will then, as in the SM, be performed by means of four unitary mixing matrices 
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l^L,R,T^L,R, where 



Ur = U\ 



(2.11) 



dn — T^fidn , 



(2.12) 



and the CKM matrix is given by 



VcKM 




(2.13) 



As argued in [6,45], the mixing of fermion zero modes with their heavy KK resonances 
induces flavour violating couphngs of the Higgs boson^, eventually leading to additional 
tree level contributions to AF = 2 processes. However it can straightforwardly be seen 
(see Appendix C for details) that these contributions are strongly suppressed in the 
model in question. Therefore we will not consider them any further. 

2.7 Impact of Higher KK Fermion Modes 

Even at tree level higher KK fermion modes affect flavour observables through mixing 
with SM fermions. Depending on the particular structure of the Yukawa interactions, 
like-charged fermions of any KK level mix with each other. The relevant three-by- 
three subsets of the infinite dimensional full rotation matrices necessarily deviate from 
unitarity. Besides that, the small admixture of higher KK fermion modes to SM fermions 
modifies their gauge couplings since SM fermions and KK fermion modes couple in 
general differently to the various gauge boson modes. This is true not only for the heavy 
KK gauge bosons, but in particular also for the Z boson, as fermions with different weak 
isospin mix with each other. 

As the KK fermion mixing appears as a consequence of EWSB, the corrections to 
the zero mode gauge couplings can be estimated to be of order 0(w^/M^j^). While 
this is a sub-leading and therefore small effect in the case of KK gauge boson couplings, 
whose flavour violating couplings are 0{1), the situation is a priori different in the case of 
flavour violating Z couplings that appear first at the 0{v'^ /M^^) level^. However we have 
checked numerically by including the first fermionic KK excitations that the KK fermion 
contribution to the Z coupling is generally suppressed with respect to the {Z^^\Z^^) 
contribution. Therefore their impact is subleading not only in the case of AF = 2 
observables, where the Z contributions are of higher order anyway, but also in the case 
of AF = 1 rare decays studied in [56]. In this context we underline that the custodial 
protection mechanism for the Zd^^d''^^ coupling discussed in Section 4.4 is effective not 

^We would like to thank Uli Haisch for bringing this issue to our attention. 
^We would like to thank Csaba Csaki for bringing this issue to our attention. 
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only for the {Z^^\ ^x'') contributions, but also for the KK fermion contributions, as the 
fermions in the model considered are embedded in P^^j-symmetric representations. 

The explicit analytic formulae for including these effects to arbitrary order in KK 
excitations clearly are beyond the scope of this paper and will be presented elsewhere. 
In the course of the present analysis we checked numerically that the effects of the first 
KK excitations on AF = 2 observables amount to at most 10% for a wide majority of 
generated data points. For completeness sake we mention that for a very small subset 
of data points, that display a large fine-tuning in some observable, the effect on that 
particular observable can be C(l). This is clearly due to the accidental suppression of 
the leading contribution which makes the sub-leading corrections relatively larger. Since 
the aim of the present work is to look for areas in parameter space with acceptable 
fine-tuning only, ignoring these points is justifiable in the framework of our analysis. 

Having at hand these results, we like to underline that the inclusion of the first 
fermion KK modes neither does modify the overall picture of AF = 2 observables in 
the model under consideration, nor does it affect the results stated in the following 
sections. In the light of other theoretical uncertainties, e. g. from the exchange of higher 
KK gauge excitations, that is also at the 10% level, or the uncertainties stemming from 
brane kinetic terms, we believe that our treatment of KK fermion modes is fully sufficient 
for the time being. 

3 Geometric Origin of Masses and Mixing Angles 

3.1 Analogy between RS and Froggatt-Nielsen Scenarios 

The aim to explain the observed hierarchies in the fermion masses and their flavour 
mixing matrices by making use of approximate flavour symmetries traces back to the 
late 1970s, to the well-known work of Froggatt and Nielsen [35]. In that pioneering paper 
a global U{1)f flavour symmetry has been introduced, under which the various quark 
fields carry different charges while the SM Higgs H is neutral under U{1)f. In order to 
allow for non-vanishing flavour mixing, the flavour symmetry is spontaneously broken 
by the VEV of a scalar $, the so-called flavon field, that transforms as gauge singlet, 
but is (singly) charged under U{1)f- In order to obtain small flavour violating effects 
consistent with observation, the flavon VEV ($) has to be much smaller than its mass 
m$ ~ A. The effective flavour violating parameter is then given by e = ($)/A ^ 1. 

A very similar structure is encountered in RS models with fermions living in the 5D 
bulk. In that case the flavour U{1)f symmetry corresponds to translations along the 
extra- dimensional coordinate < y < L, under which the metric is self-similar. The 
Higgs field, living on the IR brane, is external to this self-similarity of the bulk. The 
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fermions, on the other hand, are locahsed along the extra dimension by means of their 
bulk mass parameters CQ^u,di i- e. the bulk mass parameters can be interpreted as charges 
under self-similarity transformations. Self-similarity is broken explicitly by the presence 
of the IR brane, giving rise to the symmetry breaking parameter e~^^ <^ 1. The one-to- 
one correspondence between a Froggatt-Nielsen flavour symmetry and bulk fermions in 
RS is summarised in Table 1. 



Proggatt-Nielsen symmetry 


bulk fermions in RS 


CFT interpretation 


U{1)f symmetry 
U{1)f charges Qp = ai,bi,di 

VEV of scalar ^ {Qp = I) 
e= ($)/A < 1 


self-similarity along y 
bulk mass parameters Cq 
IR brane at y = L 
warp factor e~^^ 


scale invariance 

anomalous dimensions ji 
of fermionic operators 

CFT condensate 
scale ratio Ajr/Mpi 



Table 1: Correspondence between Froggatt-Nielsen flavour symmetry, bulk fermions in 
RS and their dual CFT description. 

It is also interesting to consider this correspondence in the CFT holographic dual. 
In that picture the 5D bulk coordinate y corresponds to the energy scale of the CFT. 
Translational invariance along the bulk then translates directly into scale invariance. 
This scale invariance is spontaneously broken by a CFT condensate at the Air ~ 1 TeV 
scale, corresponding to the IR brane at y = L. The symmetry breaking parameter is 
therefore given by Air/Mpi <^ 1, where the large hierarchy between these two scales 
arises naturally as the spontaneous breaking of the CFT is due to radiative corrections. 

As the Higgs field is part of the conformal sector, it couples strongly to compos- 
ite fermionic operators of anomalous dimension 7j, one for each quark species. Those 
fermionic operators then mix with elementary fermions that are external to the CFT 
and correspond to the SM quarks. The size of this mixing is directly related to the 
effective Yukawa couplings and depends exponentially on the anomalous dimensions 7^. 
Thus effectively the 7j can be interpreted as different flavour "charges" . 

The structure of the effective Yukawa coupling matrices in RS then turns out to be 
completely analogous to that analysed by Froggatt and Nielsen [35]"^. This analogous 
structure can now be used to derive analytic expressions for the quark masses and flavour 
mixing matrices Ul.r, T^l,r in terms of the fundamental model parameters. Therefore 
we have checked the respective analytic expressions in [35] and carefully adapted them 
to the RS scenario in question. The result is summarised in the next section. 

•^This has also been noticed and worked out independently in [45]. 
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3.2 Quark Masses and Flavour Mixing Made Explicit 

Assuming an IR brane localised Higgs boson the effective 4D Yukawa couplings Y'^''^ 
can be written in terms of the fundamental 5D Yukawa couplings A"'*^ and the fermion 
shape functions f^, ff {i = 1, 2, 3) as given in (2.10), where the hierarchies in the 4D 
Yukawas arises through the hierarchies 

/l'^«/2'^«/3^, (3.1) 

/r«/2"«/3", (3.2) 

/'«/2«/3- (3.3) 

Keeping only the leading terms in the hierarchies jQ''^^'^ j jQ''^^'^ <; obtain for 

the quark masses^ 



rrib 



rrid 



kL 



Q fd 



II \d- \d _ \d \ 
^33^22 ^23^32 



^33 



72 

V det(A'= 



/2^/2' : 



-^33-^22 ~ ■^23-^32 



kL 

kL 



(3.4) 

(3.5) 
(3.6) 



and analogous expressions for the up-type quark masses mt^c,u, with only replacing "A*^" 
by "A"" and 7""' by 

Similarly, for the flavour mixing matrices Ul,RiT^l,r defined in (2.11), (2.12) we find^ 



L)ij 



and 



R)ij 





(« < J 


< 1 


{i = 3 




(^ > J 


„d fi_ 
Hij fd 


(^ < 3) 


< 1 


{i = 3 


fd 


(^ > j) 



(3.7) 



(3.8) 



Analogous expressions hold for Ui ji with replacing "c?" by 



"'We would like to thank Katrin Gemmlcr for checking all formulae given in this section and Ap- 
pendix A. 

^See Appendix A for the explicit formulae of ujfj and pfj. 
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Finally, making use of Vckm = I^I^^Li "we obtain 

jQ jQ jQ 

Vus = «12-^ , Vub = "13-^ ' ^cb = "23-^ , (3.9) 

A A A 

with 

a., = J2i<y<- (3-10) 

k=i 

We would like to stress that the formulae given above are valid at leading order in 
^Q,u4 1 jQ,u4 ^ ^--j^ 1^^^ exact in the entries of the 5D Yukawa couplings A"''^. 

Finally, we comment on the complex phases in the above formulae. It can straight- 
forwardly be seen that in general the quark masses, as given in (3.4)-(3.6), are complex 
quantities. In order to obtain positive and real values for the quark masses, the unphys- 
ical phases in (3.4)-(3.6) have to be removed by suitable phase redefinitions, which will 
then also affect the phases of the flavour mixing matrices Ul,Ri T^l,r- Similarly, suit- 
able phase redefinitions have to be performed in order to work with the standard phase 
convention for the CKM matrix [46]. 



3.3 Discussion 

The formulae derived above considerably improve the widely used (see however [45]) 
naive estimates 

<''~^A^/^/r', (3.11) 

where A is the average value of the (anarchic) 5D Yukawa couplings, and 

fQ £U £d 

(Wi).„ (Di)., ~ ^ , {Unh-^, (I^Rh-jl (^<J)- (3.12) 

These estimates are obtained from assuming a completely anarchic, i. e. structureless, 
Yukawa coupling matrix. However, a random 3x3 complex matrix generically does not 
have all entries of equal size, unless there is some symmetry enforcing such a structure. 
In addition the formulae (3.11) and (3.12) give no hint how the complex phases present 
in the fundamental Yukawa couplings A"''^ are related to the effective CP-violating phases 
of the Ul,Ri 'Dl^h flavour mixing matrices and to the KM phase. 

In eq. (3.4)-(3.9), on the other hand, the dependence of the quark masses and mixing 
matrices on the elements of the fundamental Yukawa couplings A"''^ is spelled out in 
explicit terms, and the only approximation made is the neglect of higher order corrections 
in the hierarchies jQ'^^'^ ^ jQ'^^'^ (^^ < jn general, this approximation is well justified, 
as the fermionic shape functions exhibit a strong hierarchy, with the weakest one in the 
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right-handed down sector. Therefore the largest uncertainties of the above formulae are 
generally to be expected in Dr. In addition, as the exact dependence on A"''^ is calculated, 
predictions can be made not only for the absolute size of the IAl,r,'Dl^r elements, but 
also for their complex phases, that are relevant for CP-violation. 

However it can happen that the leading order terms given above are accidentally 
suppressed by the structure of the Yukawa couplings A"'*^. In such a case, the naively 
expected accuracy of the above formulae is lost, and higher order corrections have to be 
included in order to obtain meaningful results. 

In addition, the unitarity of the flavour mixing matrices Ul,Ri 1^l,r, and therefore also 
of the CKM matrix, is intrinsically violated by contributions suppressed by f^''^''^/ jQ'"''^ 
{i < j). The exact unitarity of these matrices, however, is very important for the study 
of FCNC processes. 

Therefore the above formulae should not be used to perform exact calculations of 
flavour violating observables, but are meant to give an improved estimate of the size of 
effects to be expected. 

In our numerical analysis we will make use of the formulae (3.4)-(3.9) in order to fit 
the SM quark masses and CKM mixings, which will subsequently be checked numerically. 
For the study of the AF = 2 observables in question we will then use the exact numerical 
results for the V^ ji mixing matrices. 



4 AF = 2 Transitions 
4.1 Preliminaries 

In what follows we will use conventions and notation of [47] so that an easy comparison 
with the SM results and with the results obtained in the Littlest Higgs model with 
T-parity (LHT) will be possible. 

The SM Hamiltonians for — and i?^^ — mixings, in the notation also 
used in the present paper, can be found in (3.1) and (3.2) of [47], respectively. The 
SM contribution to the off-diagonal element M12 in the neutral i^Tand 3^ meson mass 
matrices is then given as follows 

(^i'^)sM = ^^FIEktukMI, [\f^,S, + Xf^^St + 2A:A*r/3^ct] , (4.1) 



127r2 



(4.2) 



where Aj = V*yid and X^f* = Vf^Vtq with Vij being the elements of the CKM matrix. 
Here, Sc, St and 5*^ are the one-loop box functions for which explicit expressions are 
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Figure 1: Tree level contribution of KK gluons to — mixing. 

given e.g. in [47]. The factors rji are QCD corrections evaluated at the NLO level 
in [48-52]. Finally Bk and Bb^ are the well-known non-perturbative factors. The 
amplitude (M^2)sm ^^'^ obtained from (4.2) by simply replacing d by s. 
It should be emphasised that in the SM only a single operator 

{sd)v-A{sd)v-A = [s7m(1 - l5)d\ ® [si^il - l5)d] (4.3) 

and 

{bq)v-A{bq)v-A = [h,{l - ® [h^{l - 75)?] (4.4) 

contributes to and Mfg {q = d,s), respectively. Moreover complex phases are only 
present in the CKM factors. 

Our next goal is to generalise these formulae to include the new tree level contribu- 
tions from KK gluons as shown in Fig. 1. We will see that three distinct new features 
will characterise these new contributions: 

1. The flavour structure will differ from the CKM one. 

2. FCNC transitions will appear already at the tree level as opposed to the one-loop 
SM contributions in (4.1) and (4.2). 

3. In addition to {sd)v~A{sd)v-A and (bq)v-A(bq)v-A (with q = d,s) new operators 
will be present in the effective Hamiltonians in question. 

We recall that only the first feature is present in the LHT model. 
4.2 Tree Level KK Gluon Contributions 

We begin our discussion with the tree level exchanges of the lightest KK gluons G^^'*" as 
shown in the case of AS* = 2 transitions in Fig. 1. Analogous diagrams contribute to 
B^g — B^^ mixings. We will analyse tree level EW contributions subsequently. 
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The neutral current Lagrangian describing the interaction of the hghtest KK gluons 
G'^(^) (a = 1,...,8) with SM down quarks {i = 1,2,3) given first in the weak eigenstate 
basis is flavour diagonal and given by 

-CSc" = -Puv Yl ^i^/^^" l^d^PL + eR{i)PB] i'iG^a^''> , (4.5) 

i 

where puv parameterises the influence of brane kinetic terms on the SU (3)c coupling, 
as introduced in (2.8). The colour matrices satisfy [ta,4] = ifabctc- We suppress the 
quarks' colour indices (a,/3) for the moment, Pr^l = (1 ± 75)/2- The value of puv is 
very important for the present analysis and is model dependent as discussed in [20]. We 
recalled this discussion in Section 2.5. 

Siii) and £i^(^) are given by the overlap integrals of quark shape functions and the 
shape function of the G'^^^\ with the latter strongly peaked towards the IR brane. From 
(2.8) we have 



1 



L 



eL,R{i)=9T^j^ dye'y[ff^^{y,d^)\ g{y) , (4.6) 

with f^^jiiy, c\,) given in (2.6), (2.7) and g{y) in (2.3). As the shape functions of L and R 
quarks generally differ from each other, £l(z) 7^ Bfi{i) and parity is broken by QCD-like 
interactions in this model. Most importantly eL,R{i) depend on the flavour index i. This 
dependence breaks the flavour universality of strong interactions and implies tree level 
FCNC transitions mediated by G^*-^^ as we will see in a moment. 

In what follows it will be useful to introduce two diagonal matrices 

Sl = diag(£i(l),£z.(2),£L(3)) (4.7) 

= diag(£R(l),£R(2),eR(3)) (4.8) 

and two non-diagonal matrices 

AL = V[eLVL (4.9) 

An = V^j,inVn (4.10) 

with Vi^fi defined in (2.12). A^^r describe the FCNC couplings of down quark mass 
eigenstates to the lightest KK gluons. 

After rotation to the mass eigenbasis we find then 



^QCD ^ (^^Tm^'^^l) + Al^hi.f^dL) + A'ahl.t^SL)] Gf^ , (4.12) 
^QCD ^ \^^sd^sR^,t'^dn) + A'^ihn^.edn) + A^ibni^f^sn)] G^^^'^ , (4.13) 



.QCD . ^QCD 



where 
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and A|^^ are the elements of the matrices A^ /j. These elements are complex quantities 
and introduce new flavour and CP-violating interactions that can have a pattern very 
different from the CKM one. 

The diagrams in Fig. 1 lead to the following effective Hamiltonian for AS = 2 tran- 
sitions mediated by the lightest KK gluons with mass Mkk: 

+ 2At^t {sLl.t'^dL) {sRl^edR) ] . (4.14) 

For the B^^ — B^^ Hamiltonians one has to replace "s(i" by "M" and "6s", respectively. 

The Hamiltonian in (4.14) is valid at scales O(Mkk) and has to be evolved to low 
energy scales fi = 0(2 GeV), ^{rrif,) at which the hadronic matrix elements of the opera- 
tors in question can be evaluated by lattice methods. The relevant anomalous dimension 
matrices necessary for this renormalisation group evolution have been calculated at two- 
loop level in [53, 54] and analytic formulae for the relevant QCD factors analogous to 
rji in (4.1) and (4.2) can be found in [55]. When using these formulae we neglect the 
unknown 0{as) corrections to the Wilson coefficients of the relevant new operators at 
H = Mkk- As ^^(Mkk) with Mkk — 3TeV is small, these corrections are negligible in 
comparison with the effects of renormalisation group running from /i = Mkk down to 
/i ~ 0(m6,2GeV). 

Our next task is then to transform the operator basis in (4.14) into the basis used 
in [55]: 

Qr'' = is^f^Pid) (s^PLd) , 
= (si.PRd) (s^PRd) , 
Qf« = (s^^PLd) (sYPRd) , 

Q}f = {sPLd) {sPRd) , (4.15) 

where we suppressed colour indices as they are summed up in each factor. For instance 
si^Phd stands for Sal^Pida and similarly for other factors. 

A straightforward calculation gives us the effective Hamiltonian^ for AS* = 2 transi- 
tions in the basis (4.15) with the Wilson coefficients corresponding to fi = (9(Mkk) 

[njAS=2-\ ^ Ir^VLLriVLL I r^VRRriVRR 

[Hes Jkk = 7772- 1^1 ^1 ^1 
^J^KK 



^C^RqLR + c^RqLR^ ^ (4.16) 



^Wc would like to thank Michaela Albrccht for checking this resuh. 
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where 



C^r^(MKK) = ^Puv^ (AS' 

C^\Mkk) = -Ap^v'^AtAt . (4.17) 

Analogous expressions exist for — and B^ — B^ systems, with sd replaced by bd 
and bs, respectively. We confirm the results of [20]. 



4.3 Tree Level Electroweak Contributions to AF = 2 Processes 

The KK gluon tree level contributions in Fig. 1 discussed until now are believed to 
dominate the NP contributions to AF = 2 processes in the model in question. However 
we will demonstrate now that while this is justified in the case of Ek and AM^, in the 
case of Bd^s physics observables it is mandatory to include also tree level EW gauge boson 
contributions. To our knowledge the only paper studying EW contributions to AF = 2 
processes is the analysis of Burdman [31]. However, in that paper only Z contributions 
have been considered, and moreover the QCD renormalisation group enhancement of the 
Qlr operators has not been taken into account. As will be demonstrated below in the 
model considered here the dominant EW contributions do not come from Z but from 
tree level exchanges of the two new heavy gauge bosons Zh and Z'. The contribution of 
the KK photon turns out to be much smaller than the latter contributions. 

Let us begin with the KK photon contribution A^^\ The contributing diagrams are 
as m Fig. 1 with replaced by A^^\ is given simply by (4.5) with the colour 

matrices t° replaced by the identity in colour space and £:l,_r(0 replaced by the overlap 
integrals similar to the ones in (4.6) but with gf^ replaced by e^°, the 4D QED coupling 
constant. Note that the shape function of the KK photon is equal to g{y)- Moreover it 
is useful to absorb the electric charge factor in eL^R{i){A^^^). 

The calculation is simplified relative to the previous one by the fact that in the 
absence of one immediately obtains the result in the basis (4.15). 
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We find the following corrections to the Wilson coefficients Cj(MKK) 

2[Af(A«)]% 

0, (4.18) 

where A^£'^j^{A^^'') are analogous to A^£'^j^{G^^^) considered before with explicit expressions 
given in Appendix B. The following observations should be made: 

• A'^fj^iA^^^) are suppressed by the charge factor 1/9 and «QED/tts(^-^KK) relatively 
to A2'^(G^^^). These suppressions are partially compensated by the absence of the 
1/3 colour factors in (4.18). 

• Without 0{as) corrections to the tree level exchange of the KK photon, the coef- 
ficient [AC|''^(Mkk)] "^^^ vanishes. Strictly speaking for a NLO-QCD analysis the 
0{as) corrections to the result (4.18) should be included. But as these corrections 
are small we can neglect them. 

• The mixing of with Qf^ generates through renormalisation group effects 
a non- vanishing [AC|'^(/io)] ^™ that is proportional to [ACi^{Mkk)]^^^ and 
consequently is C(aQED)- 

We consider next the contributions of Z, and Z' gauge bosons that before EWSB 
correspond to the zero mode Z^^\ its first excited KK state Z^^-* and the heavy 
gauge boson, a linear combination of W'^ and of U{l)x [37]. Clearly while the KK 
gluon and photon contributions are universal to all RS models with bulk fermions, the 
contributions discussed in the following depend sensitively on the EW gauge group and 
the choice of fermion representations. 

Before EWSB the couplings of Z^^^ to quark flavour eigenstates are universal but 
the ones of Z^^^ and Z^^ are not. After EWSB the mixing between Z^^^ and 

implies breakdown of universality of the couplings of the mass eigenstates Z, Zh 
and Z' to quark flavour eigenstates which after the rotation to quark mass eigenstates 
implies tree level FCNC processes mediated by these three gauge bosons. Now the FCNC 
couplings of Zh and Z' are 0{1), while the ones of Z are 0{v'^/M^^). Consequently, 
its contribution to AF = 2 processes is 0{v^/M^^) and can be safely neglected already 
for this reason. In addition as we will demonstrate in Section 4.4 the flavour violating Z 
coupling to left-handed down-type quarks vanishes in the limit of exact symmetry. 



[ACr^^(MKK)]^"" = 
[ACr^^(MKK)]'''^° = 
[ACi^^(Mkk)]^''° = 
[AC,^^(Mkk)]^'^° = 
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so that these contributions are suppressed in the model in question also in the case of 
AF = 1 processes [56]. 

The calculation of (!?(f^/M^j^) tree level contributions from Zh and Z' proceeds 
similarly to the calculation of the A^^^ contribution and we find: 



EW 



[AC1'^^(Mkk 



Af(Z«)) +(Af(zj;))) 



EW 



[ACf^(MKK 



(4.19) 



where the overlap integrals A^^'^^i^Z^^^) and A^''^(zj^'') are explicitly given in Appendix 
B. They include the relevant weak couplings and weak charges. 

In order to estimate the size of EW contributions when compared to the KK gluon 
exchanges we factor out all the couplings and charge factors from A|^^. The remaining 
^^LR then universal for all the gauge bosons considered up to the different boundary 
condition of zj^^ on the UV brane, whose inclusion amounts to only a percent effect on 

Adding the contributions (4.17), (4.18) and (4.19) and evaluating the various cou- 
plings we find 



Cr^(MKK) 



sd\2 



(0.67 + 0.02 + 0.56)(Af )^ = 1.25(A^ ; 
(0.67 + 0.02 + 0.98)(A^"')2 = 1.67(A^'^)2 
(-0.67 + 0.04 + 1.13)(Af A^^) = 0.50(Af A^'^) 



(4.20) 



where the three contributions correspond to KK gluon, KK photon and combined {Z', Zh) 
exchanges respectively.^ The Wilson coefficient C|^^(Mkk) receives only KK gluon con- 
tributions at /i = Mkk- 

We observe that the EW contributions are dominated by Z', Zh exchanges and in 
the case of CY^^, CY^^ and Cf^ amount to +87%, +150% and -175% corrections. In 
particular the sign of Cf^(MKK) is reversed. 

We conclude that the EW gauge boson contributions to the Wilson coefficients 
CY^^ , CY^^ and Cf ^ at yU = Mkk are of the same order as the KK gluon contributions 



^These results are obtained neglecting the running of the EW gauge couplings between the EW scale 
Mz and the KK scale Mkk- Taking into account also these contributions, we would have corrections 
to the gauge couplings at the 5% level, so that we can easily neglect them. 
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and have to be taken into account. In the case of Ek and AM^ the strong enhance- 
ment of the coefficient C|'^ through QCD renormahsation group effects and the chiral 
enhancement of the hadronic matrix element of Q2 ^ assure that KK gluon contributions 
still dominate by far over EW contributions, although the reversal of the sign of Cf^ 
makes the constraints from Ek and AM^ to be slightly stronger. 

However, in the case of B^^s physics observables the QCD renormahsation group 
enhancement in the LR sector is smaller than in the K sector and the chiral enhancement 
of (Q2^) and (Qf^) is absent. Therefore the Q!(^^ operator becomes important even 
without the EW contributions and it is further enhanced when these contributions are 
taken into account. 

At first sight our finding that EW contributions can compete with QCD contributions 
is surprising. On the other hand one should remember that KK gluon contributions sim- 
ilarly to EW contributions are suppressed by their large masses and the main difference 
between these contributions results from gauge couplings, colour factors, weak charges 
and renormahsation group effects. Our analysis shows that with the exception of all 
these effects conspire to make EW heavy gauge boson contributions to be as important 
as the KK gluon contributions in Bd,s physics AF = 2 observables. 



4.4 Custodial Protection of Tree Level Z Contributions 



Applying the same method already used to compute the contributions of the EW bosons 
Zh and Z' to AF = 2 processes, we find the corrections to the Wilson coefficients due 
to the Z exchange 



[ACf^^(MKK 

[ACr^^(MKK 
[AC^'^(Mkk 

[ACi^«(MKK)]^ = 0^ 



Mkk 

Mz 
Mkk 

Mz 
Mkk 



{It 



Af -rAf (Z^'^ 



X 



Af - rA|j^(z(') 



X 



{Xtf Al'iZ^'^)-rAtiZ^i^^ 



A^^(Z(^)) - rA'j^iZ^^^ 



X 



(4.21) 



where we have defined the quantity r = cos cos (p = f cos ^ cos (j), and are the 
overlaps of the gauge boson shape functions with the Higgs profile as defined in [37, 56] . 

We see explicitly that the Z contributions to AF = 2 processes are suppressed first 
by (M^/Mkk)^ relative to Zh and Z' contributions. But in fact in the case of CY^^ 
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and Cf the suppression is much stronger as the custodial symmetry relevant for the 
Zb]J)L protection is also active here, and it is violated only by the boundary conditions 
on the UV brane. Neglecting this breakdown as in the estimates of (4.20) and using the 
couplings in Appendix B [37] we find 

ACf^(MKK) ~ (1 - r)(Af )(A^'^) . (4.22) 

In the limit of exact Plji symmetry f = 1, so that these two contributions vanish. As the 
right-handed couplings of Z are smaller anyway its contribution to AF = 2 observables 
is negligible. 

We recall from the discussion in Section 2.7 that KK fermion contributions to the Z 
couplings appear at the same order in the v'^/M^^ expansion as the KK gauge contribu- 
tions already discussed and therefore have to be considered as well. It is now important 
to note that provided all quark representations in the model are symmetric under the Plr 
exchange symmetry, as is indeed the case in the model considered, the custodial protec- 
tion mechanism is effective also for the latter contribution. Again, small non-vanishing 
contributions appear due to the symmetry breaking by UV boundary conditions, but 
they are found numerically small as expected. 

These findings have also implications for AF = 1 processes. Also there the left- 
handed couplings of Z to quarks are strongly suppressed by the custodial symmetry so 
that new physics contributions to meson decays with leptons in the final state turn out to 
be dominated by tree level right-handed couplings Zd^^d?^ [56] . This should be contrasted 
with the model considered in [45] where the protection of the Zd^d'^ couplings is absent 
and tree level Z contributions to AF = 1 processes are expected to be significantly 
larger. 

4.5 Afi2 from KK Gauge Bosons 

Denoting the contributions of KK gluons to the Wilson coefficients in (4.17) by [Cj(MKK)]'^, 
we finally have 

^(Mkk) = mM^K)f + [AQ(Mkk)]'^''° + [Aa(MKK)]'''^ , (4.23) 

with the various contributions given in (4.17), (4.18) and (4.19), respectively. 

The renormalisation group evolution from /i = Mkk to a low energy scale /iq can 
be done separately from the additive SM contribution, even if Q)[^^ is equal up to a 
factor of 1/4 to the SM operator {sd)y_^j^{sd)y_j^. We recall that QX^^ and QX^^ 
renormalise without mixing with other operators and that their evolution is the same as 
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QCD is insensitive to the sign of 75. But as CY^^{Mkk) ^ Ci^^^(Mkk), their Wilson 
coefficients at /iq will differ from each other. On the other hand Q^^ and ^ mix under 
renormalisation so that the RG evolution operator is a 2 x 2 matrix. 

The outcome of this analysis is an effective Hamiltonian relevant at the low energy 
scale /Iq 

+ ^(/io) ^ + C2^^(/io) Q^^] , (4.24) 

with analogous expressions for the AB = 2 Hamiltonians. 

The contribution of the KK gauge bosons G^^\ A^^\ Zh, Z' to the off-diagonal element 
is then obtained from 

(M,^);k = [^c^/^'Ikk • (4-25) 
To this end one has to evaluate the hadronic matrix elements 

= (4.26) 

They can be parameterised as follows 

(QV'if^)) = (Qr^(/i)) = IralFlcBr'if^) , (4.27) 

(Qf^(/i)) = -^-R{^,)mlF'^B^^{^^) , (4.28) 

(Q2^(/i)) = ^^(/iXi^|52^'^(/^) , (4.29) 
where the Bi parameters are known from lattice calculations. They are related to the 
parameters Bi, B^ and i?4 calculated in [57,58] as follows 

BY'^^fi) = 5i , fif ^(/i) = B, , B^^'ifi) = B, , (4.30) 

and their numerical values are given in Table 2. It should be stressed that Bi{fi) are not 
renormalisation group invariant parameters in contrast to Bk in (4.1) but in view of the 
results in [55, 57, 58] it is easier to use them in this way. Finally 

Rif^)=( . . X ■ (4.31) 

Collecting all these results we find (/il = 2 GeV) 

~R{^^L)C^^^{^^L)B^ + ^Ri^^L)C^^{^^L)B^] * . (4.32) 
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Analogous expressions can be derived for (M{^2)kk (^i2)kk 
and — mixings, respectively. For instance (/ib = 4.6 GeV) 

1 



relevant for 



B', 



KK 



with 



(4.33) 



(4.34) 



The values of the Wilson coefficients Cj in (4.33) differ from those in (4.32) as different 
A*-' are involved and the scales fii and fib in (4.32) and (4.33) are different from each 
other. Similarly Bf in (4.33) differ from the ones in (4.32) as now hadronic matrix 



elements between B^ and B^ are evaluated. 



The values for Bi in the MS-NDR scheme that we will use in our analysis have been 
extracted from [57] and [58] for the K'^ — K^ system and — system, respectively. 
They are collected in Table 2, together with the relevant values of fiQ. 





Bi 


^4 


B, 






0.57 


0.81 


0.56 


2.0 GeV 


B'^-B'^ 


0.87 


1.15 


1.73 


4.6 GeV 



Table 2: Values of the parameters Bi in the MS-NDR scheme obtained in [57] (K^-K^) 
and [58] (B^-B^). The scale fiQ at which Ci are evaluated is given in the last column. 
For Bk m (4.I) we use Bk = 0.75 ± 0.07 [59]. 



4.6 Combining SM and KK Gauge Boson Contributions 

The final results for M^^, Mf^ and M^^, that govern the analysis of AF = 2 transitions 
in the RS model in question, are then given by 

K2 = (^12) SM + (Ml,) (z = K, d, s) , (4.35) 

with (M^2)sM given in (4.1)-(4.2) and (M*2)kk in (4.32) and (4.33). 

4.7 Basic Formulae for AF = 2 Observables 

We collect here the formulae that we used in our numerical analysis. We would like to 
emphasise that, although physical observables are phase convention independent, some 
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of the formulae collected in this section depend on the phase convention chosen for the 
CKM matrix and yield correct results only if the standard phase convention [46] is used 
consistently. 

The Kl — Ks mass difference is given by 

AM^ = 2 [Re (M,^)g^ + Re {M^,)^^] (4.36) 

and the CP-violating parameter ek by 

= ^^^M^^^^^ (^1^^) SM + (M^^) kk] ' (4-37) 

where Lf^ = (43.51 ±0.05)° and = 0.92±0.02 [60] take into account that ip^ ^ 7r/4 and 
includes an additional effect from ImAo, the imaginary part of the 0-isospin amplitude 

in ^ TTTT. 

For the mass differences in the B'^^ — B^^ systems we have 

AM, = 2 |(Mf2)sM + (K2)kkI (Q = d,s). (4.38) 
Let us then write [61] 

Mf, = (Mf2)sM + (^i2)kk = (^i2)sM Cb^'''^^ (4.39) 

where 

«2)sM=l(^i'2)sMh'^'' /^^22°, (4.40) 
(M^,)^^ = \iM^,)sM\ e''^' , Ps ^ -1° . (4.41) 
Here the phases P and Ps are defined through 



We find then 



and 



Vtd = and = -|V^t,|e-'^» . (4.42) 

AM, = {AM,)smCb, (4.43) 

S^Ks =sin(2/5 + 2^Bj, (4.44) 
S^^ = sm{2\ps\-2^Bj, (4.45) 

with the latter two observables being the coefficients of sin(AMrft) and sin(AMst) in the 
time dependent asymmetries in B^ ^ "ipKs and B^ ipcj), respectively. Thus in the 
presence of non- vanishing ^9^^ and (pB^ these two asymmetries do not measure P and Ps 
but [P + (pBa) and {\Ps\ - V^sJ, respectively. 
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At this stage a few comments on the assumptions leading to expressions (4.44) and 
(4.45) are in order. These simple formulae follow only if there are no weak phases in the 
decay amplitudes for — > ipKs and —>■ ipcf) as is the case in the SM and also in the 
LHT model, where due to T-parity there are no new contributions to decay amplitudes 
at tree level so that these amplitudes are dominated by SM contributions. In the model 
discussed in the present paper new contributions to decay amplitudes with non- vanishing 
weak phases are present at tree level. However, these new contributions are suppressed 
by M^/M^j^ and, as they involve charged currents, they can be safely neglected with 
respect to the SM tree level contributions. Basically in the present analysis we make a 
working assumption that tree level contributions from new physics can only be important 
in processes in which SM contributions are absent at tree level as is the case for and 
discussed above. 

Now in models like the LHT model, the only operators contributing to the ampli- 
tudes and are the SM ones, that is with the {V — A) ® {V — A) structure. 
Consequently the new phases ^pb^ and ^ps^ have purely perturbative character related to 
the fundamental dynamics at short distance scales. The situation in the RS model in 
question is different. As now new operators contribute to the amplitudes the pa- 
rameters Cb^ and ipB^ in (4.39) are complicated functions of fundamental short distance 
parameters of the model and of the non-perturbative parameters B\, B^ and 5^. Thus 
the test of the RS model considered with the help of particle-antiparticle mixing and 
related CP-violation is less theoretically clean than in the case of the LHT model. On 
the other hand one should also emphasise that the main theoretical uncertainty in (4.33) 
comes from Fb^ and not from the Bi parameters. 

Finally, we give the expressions for the width differences ATg and the semileptonic 
CP-asymmetries 
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(4.46) 



the off-diagonal matrix element F^2) the absorptive part of the B^ — B^ amplitude. We 
refer to Section 3.8 of [47] for further details and just quote here [62] 





(4.48) 




(4.49) 
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Finally, we recall the existence of a correlation between A^^ and S^^j, that has been 
pointed out in [63] and which has been investigated model-independently in [64] and in 
the context of the LHT model in [47]. We will see below that such a correlation also 
exists in the model considered here. 

4.8 Summary 

In this section we have calculated the contributions of tree level KK gluon and EW 
gauge boson exchanges to the amplitudes M^, and Mfj in RS models with custodial 
protection of the Zd^d^j^ coupling. We have then given formulae for AMk-, AM^, AM^, 
S^Ks^ AFg and A%^^ in a form suitable for the study of the size of the new RS 
contribution. The numerical analysis of these observables will be presented in Section 6. 

While particle-antiparticle mixing in RS models has already been discussed in the 
literature, our analysis goes beyond these papers as we performed the full renormalisation 
group analysis, calculated both KK gluon and EW gauge boson contributions, and also 
considered more observables of interest. 

5 Strategy for the Numerical Analysis 
5.1 Flavour Parameters 

Let us begin this section by counting the flavour parameters in the quark sector, following 
[6], for completeness. 

First the 3x3 complex 5D Yukawa coupling matrices 

A", A'^ (5.1) 

contain each 9 real parameters and 9 complex phases. This is precisely the case of the 
SM. 

New flavour parameters enter through the three hermitian 3x3 bulk mass matrices 

cq , c„ , Cd, (5.2) 

which bring in additional 18 real parameters and 9 complex phases. 

In total we have thus 36 real parameters and 27 complex phases at this stage. Not 
all of these however are physical and some of them can be eliminated by the flavour 
symmetry ?7(3)^ of the 5D theory which exists in the limit of vanishing A"''^ and cq^u4- 
This flavour symmetry is identical to the one present in the SM, and as in the SM 9 real 
parameters and 17 phases can be eliminated by making use of this symmetry. Note that 
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one phase cannot be removed as it corresponds to the unbroken U{1)b baryon number 
symmetry. 

We are then left with 27 real parameters and 10 complex phases to be compared to 9 
real parameters and one complex phase in the SM. Evidently the new 18 real parameters 
and 9 phases come from the three bulk mass matrices cq, Cu and q. 

As already stated in Section 2.6, it is convenient to work in the special basis in which 
the bulk mass matrices CQ^u,d are diagonal and thus comprise only 9 real parameters. 
The remaining 18 real parameters and 10 physical phases are then collected in the 5D 
Yukawa coupling matrices A" and A'^. For our numerical analysis it will be essential 
to have an efficient parameterisation of A"'"^ in terms of only these parameters. Such a 
parameterisation will be presented in the next section. 

5.2 A Useful Parameterisation of A"'*^ 

As every complex 3x3 matrix, the 5D Yukawa matrices can always be singular value 
decomposed as 

A" = e^<^"f/iD,K , A"^ = e^'^'^UdDdV, , (5.3) 

where the Du,d are real and diagonal and the Uu,d,Vu,d G SU{3). The singular value 
decomposed representation contains redundancies which we will try to get rid off in the 
following. At this stage the right hand sides in (5.3) contain each (0, 1) + (3, 5) + (3, 0) + 
(3,5) = (9,11) parameters, corresponding to 9 real parameters and 11 phases. Two of 
those phases are of course spurious (see below) since a complex 3x3 matrix should 
be described by (9, 9) parameters. In order to find a description in terms of physical 
parameters only we will use the Euler decomposition for SU (3) matrices [65] 

U{a, a, 7, c, /?, b, 9, <p) = e'^^'^e'^^^e'^^^e'^^^e'^'^'^e'^^^e'^^'^e'^^'^ (5.4) 

where a, b, c are mixing angles and a, 7, P, 9, (j) are phases. In the basis in which CQ^d,u 
are diagonal and real we have the freedom to make the following rephasing 

^ e'^^^v.ie-'^^'^^'^QL (5.5) 
Ur e-'^^e-'^-'^^-e-'^^'^^-UR (5.6) 
dR e-'^^e-'^'-'^'^cie-'^^'t'^ddR (5.7) 

The unitary matrices U,V ma, singular value decomposition are defined up to an internal 
diagonal rephasing 

UDV = ([/e^^sA+iAsB^j^^g-UsA-UsB-j^^ _ U'DV, (5.8) 
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Using this freedom and an additional rephasing of the quark fields we find the equivalence 

A" = f/i (0, au^ , 7c7„ , cu^ , (3uu , ^Uu ,0Uu,O)D^ K , ay„ , 7y„ , cy„ , , 6y„ , 0, 0) 

, av^ , 7y„ + r, cy„ , - r, , 0, r/ V3) . (5.9) 

The entries r/^/S can be again rotated to zero due to the freedom to rephase the quark 
zero modes. Using this invariance parameterised by r allows us to choose jUu = 0- We 
can now define A" and A'^ in terms of physical parameters only 

A" = Ul{0,au^,0,cu^,(3uu,bu^,0uu,^)DuVu{av^,av^,Jv^,cv^,Pv^,bv^,0,0) , (5.10) 
A'^ = Udi0,aua,-fua^cua,pua,bua,0,0) DdVd{ava,ava,'yva,cva,Pva,bva,0,0) , (5.11) 

with Du = dia.g{yl, yl, ul) and Dd = diag(?/^, yj, yj). Altogether we find 18 real param- 
eters and 10 physical phases contained in the 5D Yukawas, as expected. 

5.3 Guideline for the Parameter Scan 

The starting point of our numerical analysis is the generation of random 5D Yukawa 
coupling matrices A"''^. This can efficiently be done by means of the parameterisation 
presented in the previous section. In our scan we take 

0<<,,<3 (z = l,2,3), (5.12) 

where the upper bound stems from the perturbativity constraint on A"''^. The mixing 
angles 

cr/„, bu^, ay„, cy„, 6y„, au^, cu^, bu^, ay^, cy^, by^ (5.13) 
and CP-violating phases 

f3u^,0uu,(^Vu,lVu,Pvu,lUa, Pud, av„ Iv,, Pv„ (5.14) 

will be varied in their physical ranges [0,7r/2] and [0,27r], respectively. 

The overall scale for the bulk mass parameters Cq ^d ^i^^ be fixed by the requirement 
that Cq lies in the range 

0.1<4<0.5, (5.15) 

allowing for consistency with EW precision data thanks to the protection of the ZhiJjL 
coupling [15, 16]. The remaining bulk mass parameters will then be fitted making use of 
the analytic formulae of Section 3. 
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Having generated a 5D parameter point, we check its consistency with the measured 
quark masses and CKM parameters by diagonahsing numerically the obtained effective 
4D Yukawa coupling matrices Yu^d- Note that in order not to depend on unphysical 
phases at this stage, we choose to fit the Jarlskog determinant [66] 

JcP = Im(K.KsCK:) (5.16) 

rather than the CKM angle 7 = aigiVub)- 

Finally, in order to be able to work with the well-known formulae summarised in 
Section 4, we remove unphysical phases by proper phase redefinitions of the quark fields, 
requiring real and positive masses and standard CKM phase conventions. 

Throughout the major part of our analysis, we will keep fixed 

/ = lTeV, (5.17) 

corresponding to 

MKK^2.45TeV. (5.18) 

Only in the last part of our analysis, where we wish to determine a generic bound on 
the KK scale by the requirement of naturalness, we will also vary /. 



6 Numerical Analysis 
6.1 Introduction 

Having at hand all the relevant formulae for AF = 2 processes in the RS model in 
question, we will investigate how much fine-tuning of parameters is necessary in order 
to obtain a satisfactory description of the existing data and whether some characteristic 
patterns of deviations from MFV can be attributed to this model. 

To this end we will use the measure of fine-tuning introduced by Barbieri and Giudice 
[67] and most commonly used in the literature. In that paper the amount of tuning 
'^BciOi, pj) in an observable Oj with respect to a parameter pj is defined as the sensitivity 
of Oi to infinitesimal variations oi pj. Explicitly, 

(6.1) 

where the normalisation factor pj/Oi appears in order not to be sensitive to the absolute 
size of Pj and Oi. The overall fine-tuning in the observable Oi is then given by 

^BciOi) = maXj=i,...,„{ABG(C'i,Pj)} , (6.2) 



ABG{Oi,Pj 



Pj dOi 



Oi dpj 
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where the index j runs over all m dimensions of parameter space. Obviously, the larger 
^BciOi), the more sensitive is the value of Oj to small variations in the parameters pj, 
i. e. the more fine-tuning is required to keep Oi stable. 

Now, the RS model in question has many parameters. Moreover AM^ and AMs suffer 
from sizable uncertainties originating dominantly in Fb^ and Fb^ that enter squared in 
AM^^s and are known only within 10% accuracy. This should be contrasted with the 
CP- violating parameter ek, where the decay constant Fk is known with 1% accuracy 
and the parameter Bk, that enters Ek linearly, should be known within 3% accuracy 
already at the end of 2008 from lattice calculations with dynamical fermions [68]. Finally 
let us recall that the CP-asymmetries S^^ and S^Ks basically free from hadronic 
uncertainties and the hadronic uncertainties in the ratio AM^/AMg amount to roughly 
4%. 

On the experimental side the data on AMd, AMs and Ek are very precise, so that 
their experimental errors can be neglected for all practical purposes, while S^Ks is known 
with an uncertainty of ±4%. AM^, while very accurately measured, is subject to poorly 
known long distance contributions and we will only require that (AMx)exp is reproduced 
within ±50%. 



A = |V;,| = 0.226(2) 




Gf = 1.16637-10-5 GeV" 


-2 


\Vub\ = 3.8(4) ■ 10-3 




Mw = 80.403(29) GeV 




iKbl =4.1(1) -10-2 


[69] 


a{Mz) = 1/127.9 




7 = 80(20)° 


sin^ 9w = 0.23122 




AMk = 0.5292(9) ■ lO'^ps^^ 




= 497.648 MeV 




l^i^l = 2.232(7) ■ 10"=^ 


[46] 


niB^ = 5279.5 MeV 




AMd = 0.507(5) ps-^ 




mB, = 5366.4 MeV 


[46] 


AM, = 17.77(12) ps-^ 




r]i = 1.32(32) 


[48] 


S^Ks = 0.681(25) 


[70] 


V3 = 0.47(5) 


[49,50] 


The = 1.30(5) GeV 




V2 = 0.57(1) 




Tht = 162.7(13) GeV 




Vb = 0.55(1) 


[51,52] 


Fk = 156(1) MeV 


[71] 


Fb^ = 245(25) MeV 




Bk = 0.75(7) 




Fb^ = 200(20) MeV 




Bb^ = 1.22(12) 




^B.-y^B. = 270(30) MeV 




Bb, = 1.22(12) 




Fb,^ Bb, = 225{25) MeV 




BbJBb, = 1.00(3) 


[59] 


e = 1.21(4) 


[59] 


as{Mz) = 0.118(2) 





Table 3: Values of the experimental and theoretical quantities used as input parameters. 
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// = z (oreV 


ji = 4.0 (oreV 


fi = 1(2 (ore V 


fi = 6 lev 




3.0(10) MeV 


2.5(8) MeV 


1.6(5) MeV 


1.4(5) MeV 




o.O^io j iVleV 


4.y(izj MeV 


o.2(o j MeV 


MeV 




110(15) MeV 


90(12) MeV 


60(8) MeV 


50(7) MeV 




1.04(8) GeV 


0.85(7) GeV 


0.55(4) GeV 


0.45(4) GeV 






4.2(1) GeV 


2.7(1) GeV 


2.2(1) GeV 








162(2) GeV 


135(2) GeV 



Table 4: Renormalised quark masses at various scales, evaluated using NLO running. 
The la uncertainties are given in brackets. 

With this pattern of uncertainties in mind, we will perform our numerical analysis in 
several steps as follows: 

Step 1 

We will require that the masses of the SM quarks are reproduced within 2a. For the 
three mixing angles of the CKM matrix represented in our analysis by \Vus\, \Vub\ and 
\Vcb\, we will require agreement within 2a. As the value of the phase 7 = 5ckm from tree 
level decays still suffers from large uncertainties, we will just require that it lies in the 
range 60° < 7 < 100°. The strategy for performing efficiently step 1 has been outlined 
in the previous section. 

Step 2 

Having constrained moderately the space of parameters in the first step we will investi- 
gate how much fine-tuning is necessary in order to reproduce the experimental value of 
Ek- Similarly, we will consider the cases of AMk, being sensitive to Re{M^), and S^Ks^ 
being the most accurately known AF = 2 observable in the B systems. 

Step 3 

At this stage we will impose the experimental constraints from AF = 2 observables. In 
order not to complicate our analysis, we will set all input parameters collected in Table 
3 to their central values and instead allow the resulting observables AMk, AM^, AMg, 
AMd/AMs, Ek and S^Ks to deviate by ±50%, ±30%, ±30%, ±20%, ±30% and ±20%, 
respectively. These uncertainties may appear rather conservative, but we do not want to 
miss any interesting effect by imposing too optimistic constraints. A similar strategy for 
the error analysis has been followed in [47, 72] in the context of the LHT model. Recently 
that analysis has been updated and extended by a more careful error analysis [73], 



31 



revealing that the simphfied error analysis of [47, 72] did not have a qualitative impact 
on the results obtained. 

Furthermore additional theoretical uncertainties enter our analysis due to the several 
approximations made. First of all we are taking into account only the first gauge KK 
modes and not the full KK towers, whose neglect amounts to an error ^ 10%, which we 
have checked numerically. However one should keep in mind that the model becomes non- 
perturbative already after the first few KK modes, so that we think the result obtained 
from the sum over the whole KK tower cannot be fully trusted. In addition we do not 
numerically include the effect of mixing of the SM quarks with their heavy KK partners, 
which turns out to be subleading and at the level of 10% (see Section 2.7 for details). 

Step 4 

Having at hand those regions of the parameter space that are consistent with all available 
constraints, we study the KK gauge boson effects on those AF = 2 observables that are 
not yet known with good accuracy. These are mostly the CP- asymmetries S*^^ and ^gL, 
but also the width difference AF^/Fs, in the Bg system. 

Step 5 

Finally we will investigate whether the results obtained in Step 4 depend significantly on 
the fine-tuning Abg{£k)- Therefore we will impose the additional constraint Abg{£k) < 
20 and redo the phenomenological analysis performed in Step 4. 

Step 6 

Last but not least, motivated by the analysis in [20], we will derive a generic lower bound 
on Mkk from ek, demanding that the average required fine-tuning to get an acceptable 
Ek value does not exceed a certain naturalness limit. In this step we will therefore also 
vary the scale /. 

Throughout our analysis we will consider density plots rather than scatter plots, as these 
offer the additional information which effects are the most likely ones. We also show the 
colour bar for each of the plots, although the absolute number of points in each counting 
bin depends of course on the number of points considered and on the bin size chosen. 
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Figure 2: left: Re(M]^)ft-;^/ Re(M]^)5M and Im(Mi^)i^i^/ Im(M]^)5M; plotted on logarith- 
mic axes, right: Iie{Ml2) kk o-nd lm.{Ml2)KK! normalised to \{M12)sm\ and plotted on 
logarithmic axes. 



6.2 Results 

6.2.1 RS Contribution to M/^''^'^ 

In order to get a feeling for the size of the RS contribution to AF = 2 observables, we 
show in Fig. 2 the complex {M^)YiK and {M12)kk planes. 

In the left panel of Fig. 2, we show Im(Mj'^)KK/ Ini(M]^)sM plotted as a function 
of Re(M^)KK/ Re(M]^)sM. We observe that while Re(M]^)KK has the tendency to be 
generically somewhat smaller, albeit still competitive, with the SM contribution, the KK 
contribution to Im(M;j'^) typically exceeds the SM by two orders of magnitude. This is 
due to the suppression of Im(M]^)sM with respect to Re(M]'|)sM by roughly a factor 100, 
and leads to the generic strong constraint from ek on the KK scale identified in [20]. 
Still, already from this figure we can deduce that there exist regions of the parameter 
space for which Im(M]'|)KK ^ Iiii(M;j'^)sM5 so that agreement with the data on Sk can 
be obtained even for a scale as low as Mkk = 2.5 TeV. 

In the right panel of Fig. 2, we show Re(M*2)KK and Im(M^2)KK; normalised to 
|(Mf2)sM|- We observe that the KK gauge boson contribution tends to be of roughly 
the same size as the SM contribution, and that contrary to the SM Re(M|2)KK and 
Im(M;^2)KK are generically of the same size, so that an 0{1) new physics phase can be 
expected. For completeness we mention that the case of M{^2 is very similar to and 
we do not show it here. 

Next, we aim to analyse the importance of the various operators induced by the KK 
gauge boson exchange. Therefore in the left panel of Fig. 3 we show the ratio of the Qlr 
and Qll operator contributions to (M/|)kk- In accordance with the analysis in [20] we 
observe that the LR contribution is by far the dominant one, while the LL contribution is 
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Figure 3: The ratio of the contribution of only Qlr and only Qn to {M^)kk (l^ft) and 
{M12)kk (right), as a function of {M{^)kkI {M{^)sm (i = K,s). 

typically below 10%. The reasons for this dominance is the chiral and QCD enhancement 
of the LR operator. The contribution of Qrr turns out to be negligibly small, which 
is due to the fact that the right-handed bulk mass parameters c\ violate the flavour 
symmetry much less strongly than the left-handed Cg ones. In addition the quark 
lives closer to the UV brane than the hi one and is therefore much less sensitive to the 
flavour violation induced by KK modes close to the IR brane. 

In the right panel of Fig. 3 then the ratio of the Qlb. and Qn operator contributions 
to (M^2)kk is shown. In that case the situation differs from the K system, due to absence 
of the chiral enhancement and the weaker renormalisation group QCD enhancement. 
Indeed we find that the Qn and Qlr turn out to be competitive in size, and in most cases 
Qll even yields the dominant contribution. We note that while the Qlr contribution is 
essentially unaffected by the EW contributions, they enhance Qll by roughly a factor 
2, so that the importance of LL contributions in Bd^s physics is increased by these 
contributions. Again the contribution from Qrr to Mfg is very small. The situation in 
the Bd system is very similar and we do not show it explicitly. 

6.2.2 Fine-Tuning in AF = 2 Observables 

While deducing already from Fig. 2 the possibility to obtain ek in accordance with the 
data, we are now interested in how natural such values are. Therefore we show in Fig. 4 
the fine-tuning in ek-, ^BGi^K)-, as a function of Ek- We observe that while for generic 
values Ek / {sK)e:K.p ~ O(IOO), the fine-tuning is typically relatively small, ABcl^^i^) ~ 20, 
the average required tuning strongly increases with decreasing ek-, so that generically for 
Ek ~ (£^E')exp a fine-tuning of the order /S.-qq{ek) ~ 700 is required, i.e. the amount of 
fine-tuning increases by roughly a factor 30 — 40 when going from the generic prediction 
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Figure 4: left: The fine-tuning ^bg{^k) plotted against Ek, normalised to its ex- 
perimental value. The blue line displays the average fine-tuning as a function of 
Ek- right: The same, hut displaying only the phenomenologically interesting region 

0.1 < \EK/{eK)exp\ < 10. 

for Ek down to values in accordance with experiment. In other words, a relative fine- 
tuning at the few percent level is on average required in order to obtain Ek ~ (£^_ft:)cxp- 
However, it can also be observed that although for smaller values of Ek large fine-tunings 
become more likely, even for SM-like Ek roughly 30% of the points lie still in the range 
with small tuning, AbgI^e") < 20, so that fully natural solutions to the ''ek problem" 
can be identified. 

Let us next consider the necessary amount of fine-tuning in other AF = 2 observables. 
As examples we show here AMi<-, being sensitive to Re(M]^), and S^Ks^ being the most 
accurately known AF = 2 observable in the B systems. 

In Fig. 5 we show the fine-tuning Abg(^^e') as a function of AM^. We observe 
that the tuning is generally smaller {A-bq{^Mk) ^ 20) than in the case of Ek, and 
that the smallest average values are obtained for AMk in accordance with experiment. 
This could already be expected from Fig. 2, where we found the KK contribution to 
Re(Mj^) to be of the same order of magnitude as the SM contribution, so that generically 

In Fig. 6 we show ^^g{S^Ks) as a function of S^Ks- Also in that case the average 
fine-tuning is smallest for S^Ks in accordance with the data. In addition the overall scale 
of ^-QQ^S^Ks) turns out to be different, so that typically ^-qq^S^i^Ks) ^ 5. 

6.2.3 Full AF = 2 Analysis and CP- Violation in - Mixing 

Having convinced ourselves that in principle it is possible to obtain agreement with the 
available AF = 2 data, we are ready to perform a simultaneous analysis of all available 
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Figure 5: The fine-tuning Abg{^Mk) plotted against AMk, normalised to its experi- 
mental value. The blue line displays the average fine-tuning as a function of AM^- 




Figure 6: The fine-tuning Abg{S^Ks) plotted against S^Ks- The blue line displays the 
average fine-tuning as a function of S^Ks- 
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Figure 7: lejt: normalised to its SM value, as a function of S^^. In addition to 

the requirement of correct quark masses and CKM mixings, also the available AF = 2 
constraints are imposed, right: The same, but in addition the condition Abg{^k) < 20 
is imposed. 

constraints. To this end we now impose all AF = 2 constraints on the RS parameter 
space, as described in Step 3. The points we show in the subsequent Figures 7 and 
8 are consistent with the AF = 2 data and thus fully realistic. In order to maintain 
naturalness of the theory, the plots in the right panels of these figures fulfil the additional 
constraint Abg{£k) < 20. 

In Fig. 7 we show the semileptonic CP-asymmetry A^j^ as a function of S^^. We 
observe that while values of these asymmetries close to the SM ones turn out to be 
most likely, being a consequence of the generic relation |(Mf2)KK| ~ I(^i2)sm| observed 
in Fig. 2, we find that the full range of new physics phases ips^ is possible, so that 
— 1 < 5*^0 < 1 compared to the SM value (5'^</,)sm ~ 0.04, and also A^^ can be enhanced 
by more than two orders of magnitude relative to its SM value. In addition we observe 
that the model-independent correlation pointed out in [63] and verified explicitly in the 
LHT model in [47] turns out to be valid as well in the RS model in question. Comparing 
the left and right panel with each other we find that the imposition of the naturalness 
constraint ABG(^i^) < 20 does not qualitatively modify the results obtained, although 
the overall number of parameter points shown in the plots of course decreases. 

Finally in Fig. 8 we show the width difference AFg/r^ as a function of S^^. We 
observe that due to the correlation between these two observables, a future more accurate 
measurement of AFs/Fs could help to exclude large values of S^^. Again, comparing 
the left and right panel with each other we find that the overall number of parameter 
points shown in the plots decreases when imposing Abg(^k) < 20, but the result is not 
qualitatively modified. 
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Figure 8: left: AFs/Fg as a function of S^^. In addition to the requirement of correct 
quark masses and CKM mixings, also the available AF = 2 constraints are imposed, 
right: The same, but in addition the condition Abg{£k) < 20 is imposed. 



6.2.4 Generic Bound on M^kk 

So far in our numerical analysis we have fixed the scale / to 1 TeV, corresponding to the 
KK gauge boson mass Mkk — 2.45 TeV, as we were mainly interested in studying the 
effects on AF = 2 observables of KK modes that lie in the reach of LHC We have found 
that while it is possible to fulfil all existing constraints, in particular the one from Ek even 
without significant fine-tuning of parameters, we observed that generically a significant 
amount of fine-tuning is required in order to keep Ek in agreement with experiment. 

Finally, motivated by the findings in [20] , we aim to derive a generic lower bound on 
the KK scale Mkk- In order to achieve this we impose the constraint that the average 
fine-tuning required to obtain acceptable values for Ek should not exceed a certain value, 
i.e. ABG(£/r)av. < 10 or 20. 

In addition to our previously performed scan over the 5D Yukawa couplings and bulk 
mass parameters, we now take also /, or equivalently Mkk, as a free parameter. Fig. 9 
shows the average required fine-tuning in Ek, obtained by taking the arithmetic mean of 
Abg(£^a') of those points that fulfil the Ek constraint within ±30%, as a function of Mkk- 
We observe that ABG(^i^) decreases roughly as 1/Mkk, as expected from the dependence 
of the KK gauge boson contributions to M12. As we have seen in Fig. 4 the average fine- 
tuning for points that lie in the generic region for Ek is around 20. Therefore, imposing 
then as naturalness constraint ABG(^i<')av. < 20 we obtain as lower bound on the KK 
scale 

Mkk ^18 TeV, (6.3) 

in rough accordance with the result of [20]. We note that the bound in the latter paper 
has been obtained by the requirement that the average value for the Wilson coefficients 
respect the model independent bounds. This gives the same result as requiring the 
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Figure 9: The average required fine-tuning in Ek as a function of the KK scale Mkk- 

generic prediction to be within the bounds. 

If we were to impose instead the more stringent constraint ABG(^i<')av. < 10, that is 
often adopted in the hterature, we would find the even stronger constraint 

MKK;^30TeV. (6.4) 



Still we would like to stress again, that although this bound can be considered as a 
naturalness constraint on the theory coming from e^, we have found regions of parameter 
space which yield Ek in rough agreement with experiment without any significant fine- 
tuning for a KK scale as low as 2.5 TeV. Note however that sub-leading contributions like 
the radiatively induced brane kinetic terms for the fermions are expected to dominate 
in case of accidental numerical cancellations of the leading terms. 

Thus a natural solution to the "ej^ problem" with KK gauge bosons in the reach of 
the LHC cannot be excluded although we expect it to be radiatively unstable. 

7 Conclusions 

In the present paper we have performed for the first time the full renormalisation group 
analysis at the NLO level of the most interesting AF = 2 observables within the SU{'i)c x 
SU (2) I X SU {2)r X U{l)x X Plr model, including both KK gluon and EW gauge boson 
contributions. The protective custodial and Plu symmetries in such models allow for 
consistency with EW precision tests for KK scales as low as Mkk — (2 — 3) TeV that 
are in the reach of the LHC. As pointed out in [20] for an anarchic structure of the 5D 
Yukawa couplings much higher KK scales in the ballpark of (10 — 20) TeV are required in 
order to satisfy the ^i^-constraint in the presence of KK gluon exchanges. Our detailed 
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analysis confirms these finding, but having at hand more accurate formulae allows for a 
quantitative estimate of the fine-tuning in the 5D Yukawa couplings required to reproduce 
the quark masses and CKM parameters and simultaneously obtain consistency of the 
model with ek and other AF = 2 observables for Mkk — (2 — 3) TeV. 
The main messages from our analysis are as follows: 

1. While generally Ek values turn out to be significantly larger than its experimental 
value, we find regions in parameter space in which the experimental value of Ek can 
be reproduced without large fine-tuning. The situation is different for the other 
AF = 2 observables, where the experimental constraints are naturally fulfilled 
without significant fine-tuning. 

2. Very interestingly the EW tree level contributions to AF = 2 observables mediated 
by new and Z' weak gauge bosons, while subleading in the case of Ek and AMx, 
turn out to be of roughly the same size as the KK gluon contributions in the case 
of Bd.s physics observables. The Z contributions are of (9(w^/M|;j^) and moreover 
further suppressed by the custodial protection of Zd^cPi^. 

3. The amount of fine tuning required to satisfy the AF = 2 constraints in — 
and Bg—B^ systems is considerably smaller than in the case of AM^ and Ek- This 
is partly due to the fact that the role of the dangerous Qlr operators in AB = 2 
transitions turns out to be significantly less important than in AS = 2 transitions, 
so that the contributions of the operators Qll and Qlr to the AB = 2 observables 
are of the same order. 

4. The contributions of KK gauge boson tree level exchanges involving new flavour and 
CP-violating interactions allow not only to satisfy all existing AF = 2 constraints 
but also to remove a number of tensions experienced lately by the SM, observed in 
particular in ek, S^Ks and S^^ [60,74-76]. 

5. Most interestingly the model allows naturally for S^p^ as high as 0.4 that is hinted at 
by the most recent CDF and D0 data [77-79] and by an order of magnitude larger 
than the SM expectation, S'^,^ ^ 0.04. The strong correlation between S^^ and 

shown in Fig. 7 implies then a spectacular departure of the latter observable 
from its tiny SM value. 

6. The effects of the mixing of the heavy KK quarks with the SM quarks turns out not 
to be very important in particular in view of many parameters present in the model. 
As the fermion representations in the model in question are rather complicated [37] 
and these fermions do not contribute at tree level to AF = 2 processes we leave a 
detailed analysis of these effects to a separate publication. 
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7. As a by-product we analysed the connection of RS models to the Froggatt-Nielsen 
scenario and provided analytic formulae for the effective flavour mixing matrices in 
terms of the fundamental 5D parameters. We also presented a new useful param- 
eterisation of the 5D Yukawa coupling matrices, taking into account only physical 
parameters. 

Our detailed analysis of rare K and B decays in the model in question is presented 
in [56]. 
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A Additional Details on Quark Masses and Flavour 
Mixing Matrices 

Explicit expressions for uof- and pf- in (3.7)-(3.10) are given as follows: 

\d \d \d \d \d \d 

, _ 1 , ,d — ^33^12 ^13^32 , _ ^13 , _ ^23 I \ T\ 

^ii-^^ ^12 - xd xd _ xd xd ' '^13 -Td-' ^2^-^!-^ y^-^) 

^22^33 ^23^32 ^^33 ^33 

21 = - «) * , < = - (^^3) * - (^23) * < , ^32 = - «) * • (A.2) 



xd xd xd xd \ * / \d \ * / \d ^ 

d 1 „d I ^33^21 ^31^23 \ „d / ^31 \ „d / ^32 



Pa - ^1 P12 ~ \ xd xd \d Xd j ' - I xd ) ' P23 - \ xd 

V ^22^33 ~ ^23^32/ V ^33 / V ^33 

(A.3) 

P^i = -{pi2)\ pi = - (P?3) * - (P23) * P21 , Pi2 = - {pis) * ■ ( A.4) 

The expressions for u^j and p^-, that enter the formulae for Ul^^, are obtained by replacing 
"d" by "u". 
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B Details on Electroweak Contributions 



In the case of A^^\ ^L,R{i) in (4.6) is replaced by {i = 1, 2, 3) 



(^)(A«)=geme 



4D 



dye"' f^L^y^d^) 9{y) 



(B.l) 



with g{y) being the gauge KK shape function of A^^\ Then the 3x3 matrices Al^ji{A^^^) 
are defined by 



(B.2) 



with iL,R{A^^^) diagonal matrices analogous to (4.7) with the diagonal elements given 
by £L,ij(2)(v4(i)) in (B.l). Formula (B.2) allows then to find Af^(A(^)), A^^^^A^^^) and 
A5f,^(yl(i)). 

In order to give the expressions for Al^^i^Z^'^^) and Al,r{Z^^^) we introduce 



ed^^'^) = gfLT I dye'y fr{y,cl) g{y) 



'L 



ky M), 



dye 



ky 



(1)^ 
X > 



SRmz^i^ 



4D ^ 
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4D 1 



dy e 



ky 



4 I dye 



/f (2/,4)l g{y) 



f^R\y,c\,) m 



(B.3) 
(B.4) 
(B.5) 
(B.6) 



with g{y) being the shape function of Z^-*, that differs from g{y) due to the different 
boundary condition on the UV brane. Further 

-,4D 



9fR 



9 



cosip 



1 — sin ih 

2 3 ^ 



9 



4D 



COS tjj \3 



sin^ ^ 



,4D 9 



4D 



AD 9 



coscp 

4D 



COS I 



sm 

2 6 



-1 + -sin^ 
3 



(B.7) 
(B.8) 
(B.9) 
(B.IO) 



Here g^^ and gjP are the SU{2)i and U{\)x gauge couplings, respectively. Moreover 
sin^ ip ^ sin^ Oy/ and sin </>, cos (p as functions of ip are given by the formulae 

1 , sin (h 



cosip 



V^l + sin^ (f) ' 



sin ijj 



(B.ll) 



sm 
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and can also be found in [37]. Note that in the gauge KK sector 

9z!l - cos?/;cos0/€t° = (B.12) 

is at the basis of the protection mechanism for flavour diagonal Zbibi and non-diagonal 
left-handed down quark couplings to Z. See Section 4.4 for details. 

^l,r{Z^^^) and Ai^/^(Z^'') are then defined in analogy to (B.2) through 

Al,r{Z^'^) = V{JlAZ^'^)'^l,r (B.13) 

with a similar expression for z'^^ . El^r are diagonal matrices with their elements given 
by (B.3HB.6). 



C Tree Level Flavour Changing Higgs Couplings 

In this appendix we estimate the size of the relevant Higgs vertices by making use of the 
mass insertion approximation describing the mixing of fermion zero modes with their 
heavy KK partners. See also [45] for an alternative derivation. 

We start by considering diagrams with one heavy-light transition on a fermion line 
(denoted by +). 




As the Higgs vertex in that case contains a Pr projector, while the heavy-light mass 
insertion comes along with a P/,, the leading contribution from the 1/Mkk part of the 
fermion propagator vanishes, and only the non- leading i> jM'^Y^ contribution survives. 
When acting on the external fermion, the additional jMyj^^ results in the strong chiral 
suppression m'llMYj^. 

Let us next consider the case of a heavy-heavy transition in addition to the heavy- 
light transition already considered. Naively one may expect that, as now the external 
fermions are of different chirality, the suppression factor / Mkk is absent , so that such 
diagrams yield the dominant contribution. However, one finds 

~ ■^jkK^ )kA^ )ii TTJi J 3 MjT ' 
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i. e. this type of contribution is highly suppressed not only by the v'^/M^^ factor coming 
from the two mass insertions, but in addition receives a double chiral suppression factor 
mf^/M^j^. The origin of this strong chiral suppression is in fact easy to see: The Higgs 
boson, being confined to the IR brane, can couple only to that chirality of a given fermion 
KK mode that obeys a Neumann BC on that brane; this chirality is necessarily the one of 
the corresponding zero mode. This implies that again only the dependent parts of the 
KK fermion propagators contribute. Evaluating then the Dirac structure of the above 
diagram one ends up with the result stated above. We note that in the case of a bulk 
Higgs boson the heavy-heavy transition would contain both fermion chiralities, so that 
the mf /Ml^ suppression would be absent in that case and this kind of diagram would in 
fact yield the dominant contribution to flavour changing Higgs couplings. Consideration 
of diagrams with KK fermions contributing simultaneously on both external lines does 
not change this conclusion. 

We conclude that in the present brane-Higgs scenario Higgs contributions to FCNC 
processes are negligible in the model in question, which we have also verified numerically. 
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